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Abstract 



In this paper the stability of the Korteweg-de Vries (KdV) equation is investigated. It is shown analytically 
and numerically that small perturbations of solutions of the KdV-equation introduce effects of dispersion, hence 
the perturbation propagates with a different velocity then the unperturbed solution. This effect is investigated 
analytically by formulating a differential equation for perturbations of solutions of the KdV-equation. This differ- 
ential equation is solved generally using an Inverse Scattering Technique (1ST) using the continuous part of the 
spectrum of the Schrodinger equation. It is shown explicitly that the perturbation consist of two parts. The first 
part represents the time-evolution of the perturbation only. The second part represents the interaction between 
the perturbation and the unperturbed solution. It is shown explicitly that singular non-dispersive solutions of the 
KdV-equation are unstable. 



1 Introduction 

Since the discovery of the 1ST, a great number of nonlinear differential equations have been solved using this technique 
(see for instance ref. Q and the references therein). Most of the scientific effort has been focused on finding soliton 
solutions for these equations. As firstly observed by Scott-Russell in 1834, solitons have the property that they 
maintain their shape over long time-scales Q . From this observation the conclusion can be drawn that solitons which 
are observed in nature are stable solutions. 

In this paper we focus on the KdV-equation. The 1ST approach, as introduced by Gardner, Greene, Kruskal and 
Miura (GGMT), uses the inverse problem of the Schrodinger equation to solve the KdV-equation. In their approach a 
solution of the KdV-equation can be regarded as a time-dependent potential of the Schrodinger equation ||. GGMT 
have presented a method to compute the time-dependence of the S-matrix that is used to solve the inverse problem of 
the Schrodinger equation. Furthermore, they have shown that soliton solutions can be constructed using the discrete 
part of the spectrum of the Schrodinger equation. 

Recently, it was discovered, (by solving the KdV-equation using the continuous part of the spectrum of the 
Schrodinger equation) that the KdV-equation has also singular solutions Q . It is also shown in ref. Q , that the time- 
evolution of these singular solutions can be associated with a positive Lyapunov exponent. This implies that singular 
solutions of the KdV-equation can exhibit unstable behavior. 

It is remarkable that although the KdV-equation is one of the better studied equations in the field of mathematical 
physics, the matter of the stability of this equation has not been investigated. The classical argument that explains 
that solitons maintain their shape over long time-scales is the presence of the nonlinear term in the KdV-equation 
that cancels the effects of dispersion. However in nature, solitons are always contaminated with noise and it is not 
clear that when noise is added to a pure soliton solution, the soliton retains its identity. This matter is related to 
the stability properties of solutions of the KdV-equation. In this paper we show that for the KdV-equation solitons 
posses stable behavior if the propagation velocity of the perturbation differs from the perturbation velocity of the 
unperturbed soliton. In the case of the singular solution this is also the case, but in contrast to the soliton case, the 
amplitude of the perturbation can grow dramatically. This implies that the shape of the singular solution changes. 

This paper has the following structure. In Sec. 2 we discuss the balance between nonlinearity and dispersion for 
non-dispersive solutions of the KdV-equation. It is shown that for small perturbations this balance is disturbed, and 
that the perturbation has a different propagation velocity than the unperturbed solution. In Sec. 3, we derive a general 
expression for a perturbation of a solution of the KdV-equation. It is shown that the time-evolution of the perturbation 
consists of two parts. The first part represents the time-evolution of the perturbation in absence of the unperturbed 
solution. The second part represents the interaction between the soliton and the perturbation. In Sec. 4, we give an 
expansion of the time-evolution of the perturbation in a series solution. This enables us to formulate an expression for 
the stability of solutions of the KdV-equation. In Sec. 5, an numerical example is given to illustrate that the results 
that are derived in previous sections are also valid in more general cases. The results are discussed in Sec. 6. Technical 
matter concerning the Gelfand-Levitan-Marchenko equations are added in a appendix. 
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2 Nonlinearity versus dispersion 



Consider the KdV-equation: 

U xxx ~ 6uu x + u t = . , 

u(x, t = 0) = uo(x) 

In Eq.(|l|), uo(x) represents the initial condition. Non-dispersive solutions of Eq.(Q) can be obtained by searching for 
solutions u(x, t) — u(x — ct) = u(z). In this special case Eq. (jl|) reduces to: 

u — 6uu — cu = (2) 

The KdV-equation ([[]) is in the special case of non-dispersive solutions reduced to a third-order ordinary nonlinear 
differential equation. In Eq.(|^), the derivative u' stands for j^u(z). We can reformulate Eq.(||) in the following way: 

d z (u" - 3u 2 -cu) =0 (3) 

By integrating Eq. (|J) , we find that Eq. (||) is equivalent with: 

u" — 3u 2 — cu + m = (4) 

In Eq.(|J), the constant m £ M is an arbitrary integration constant. We can perform a further simplification by 
multiplying both sides of Eq.(^) with a factor u'. The total result can be formulated in the following way: 

dz (Uu' f-u :i -]-cu 2 + mu\ =0 (5) 



X 2 

By performing a further integration we find that Eq.([s]) is equivalent with: 

i(u') 2 - u 3 - icw 2 + mu + ?i = (6) 

In Eq.(||), the constant n G M is another arbitrary integration constant. In the special case that m = n = 0, Eq.(|^) is 
equivalent with: 

u 1 = ±u^/2u + c (7) 

We have obtained that in the case of non-dispersive solutions the KdV-equation (||) is equivalent with Eq. ((?]). Eq.(Q) 
is an ordinary first-order nonlinear differential equation that can be solved directly. This leads to the following result: 

4c e ±z^ 

< z ) = ~ ( 8 ) 

2 ,,+zt/c 

A) 

In Eq.(||) the constant Do £ 1R can be chosen arbitrarily and has to be determined by the initial condition. If the 
constant Dq is negative, the solution (||) can be written as a soliton solution: 

u(z) = — -sech 2 (zy / c + xq) (9) 

If the constant Dq is positive, this reduction does not take place. It follows from Eq.(||) that in this case the solution 
of the KdV-equation has a singularity, but the solution propagates without dispersion. It has been argued in ref.Q 
that these singular solutions can be associated with a time-evolution having a positive Lyapunov exponent. This leads 
to the conclusion that apart from the stable soliton solutions, the KdV-equation has also unstable solutions. 

Soliton solutions are understood to be stable solutions of the KdV-equation because the effects of dispersion and 
nonlinearity cancel each other. It is tempting to choose Eq.(^) as a starting point for the stability analysis for the KdV- 
equation since Eq.(^) is an ordinary differential equation on which all the standard techniques of stability analysis can 
be applied. However, due to the fact that Eq.(Q) is only valid for non-dispersive solutions the corresponding stability 
analysis is only valid for solutions of the KdV-equation (including a perturbation) that maintain their shape. For non- 
dispersive solutions of the KdV-equation, the effects of dispersion and non-linearity are in balance. If a perturbation 
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is added to a non-dispersive solution of the KdV-equation, the balance between the dispersion and non-linearity is 
disturbed. In order to investigate this effect we rewrite the KdV-equation in the following way: 

u t — —cu x + (cu x - u xxx ) + 6uu x (10) 

If the last two terms on the right-hand side of Eq.(|l(i|) are removed, we obtain: 

u t = -~cu x (11) 

Eq.([ll]) has the non-dispersive solution u(x, t) = g(x — ct). If the KdV-equation has non-dispersive solutions, the 
second and the third term on the right-hand side of Eq.([l0|) cancel each other. The cu x — u xxx term on the right-hand 
side of Eq.(|l0|) describes the dispersion of the solution u(x,t). The 6uu x term on the right-hand side describes the 
effects of the nonlinearity. 

The solid line in Fig. la the solid line represents a soliton at time t = 0. The dashed line in Fig. la, represents the 
soliton that is contaminated with a 10% amplitude error. In Fig. lb, the balance between the effects of nonlinearity 
and the effects of dispersion is shown for the unperturbed soliton. The short-dashed line in Fig. lb represent the effect 
of dispersion as given by the cu x — u xxx term on the right-hand side of Eq. ( |Toj ) . The long-dashed line describes the 
nonlinearity as given by the 6uu x term on the right-hand side of Eq. (Jl^) . The sum of these curves are given by the 
solid line in Fig. lb. It follows from Fig. lb that the dispersion and nonlinearity are in balance. This is consistent with 
the fact that that the soliton propagates with velocity c while maintaining its shape. 

If the soliton is contaminated with a small amplitude error, the dispersion and nonlinearity are no longer in 
balance. In Fig.lc, the examples of Fig. lb are repeated for the contaminated soliton shown by the dashed line in 
Fig. la. Similarly as in Fig. lb, the short-dashed line in Fig.lc represents the effect of dispersion as given by the 
cu x — u xxx term on the right-hand side of Eq.(|l0|). The long-dashed line describes the nonlinearity as given by the 
6uu x term on the right-hand side of Eq.(|To|). The solid line in Fig.lc represents the sum of the dispersion and the 
the nonlinearity. An imbalance between the dispersion and the nonlinearity is introduced due to the amplification of 
the nonlinearity. As a result of the fact that the nonlinearity and the dispersion no longer cancel, Eq. ( |Tl| ) (which is 
only valid for non-dispersive solutions of the KdV-equation) , does not describe the time-evolution of the perturbed 
solution. The propagation effects of the contaminated soliton is given by Eq. ([To|) . The imbalance between the effects 
of dispersion and nonlinearity is given by the solid line Fig.lc, which is equal to the sum of the last two terms on 
the right-hand side of Eq. (|Io|) . In Fig. Id, the time-derivative of the perturbation is plotted. It turns out that the 
time-derivative of the perturbation is not equal to a constant times the space derivative of the perturbation. This 
implies that the time-evolution of the perturbation is dispersive. In the following sections this behavior is investigated 
analytically. Finally, in Fig.le, the finite difference solution of contaminated soliton is given after at t = 2.2 sec. It is 
observed that as a consequence of the imbalance between the nonlinearity and dispersion, the perturbation propagates 
with a different velocity than the unperturbed soliton, and during the course of time, the contaminated soliton takes 
its most natural form. 

From this experiment, we can conclude that as a result of the imbalance between the effects of nonlinearity and 
dispersion, the perturbation of a solution of the KdV-equation propagates with a different velocity. As a result of this 
the soliton and the perturbation separate during the course of time. 

In Fig.2ab the examples of Fig.l are repeated for the singular solution. The solid line in Fig. 2a describes the 
singular solution of the KdV-equation at t = 0. The short-dashed line in Fig. 2a describes the effects of dispersion 
as given by the cu x — u xxx term on the right-hand side of Eq. (|l0|) . The long-dashed line describes the effects of 
nonlinearity as described by the 6uu x term on the right-hand side of Eq. ([!(]) . Similarly as for the soliton, it turns out 
that also for the singular solution effects of dispersion and nonlinearity arc in balance. If a small amplitude error is 
made, this balance is no longer present. In Fig. 2b, the effects of dispersion and nonlinearity arc plotted. The solid 
line in Fig. 2b represents the sum of the nonlinearity and dispersion in the contaminated case. It turns out that the 
effects on nonlinearity and dispersion are no longer in balance. As a result of this, the perturbation propagates in the 
opposite direction the unperturbed wave. 

From the simple examples in this section, we can conclude that if nondispersive solutions of the KdV-equation 
are contaminated with small perturbations, dispersion effects are introduced. The imbalance between the dispersion 
and nonlinearity is in the coordinate frame that moves with the unperturbed solution very close to the x-derivative 
of the perturbation. This implies that the perturbation has a non-zero velocity in this reference frame so that the 
perturbation separates with time from the unperturbed solution. 
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3 The stability of localized solutions 



The result in the previous section indicates that perturbations on the initial condition of the KdV-equation propagate 
with a different velocity than the unperturbed solution. In this section the effects of different perturbations are 
investigated in a more general fashion. If a perturbation u(x,t) — ► u(x,t) + f(x,t) is substituted is Eq. ([!]), the 
following differential equation for the perturbation f(x,t) can be derived: 

fxxx ~ 6 (uf x + f x u + ff x ) + f t = , . 

f(x,t = 0)=f (x) {iZ) 

Eq.(^) represents a differential equation for the perturbation f(x,t), which depends on the unperturbed solution 
of the KdV-equation u(x, t). Eq. (|f2|) can be solved using an inverse scattering technique if a satisfying Lax-pair is 
constructed. However, because of the fact that the perturbed solution u(x,t) + f(x,t) satisfies the KdV-equation, the 
solution f(x,t) of Eq.@ can be computed directly using the techniques described in ref. ||. In Appendix A, a brief 
overview of these methods is given: 

As a starting point, we assume that the reflection coefficient corresponding to the unperturbed initial condition 
uq(x) undergoes a perturbation: 

R(k, t = 0) -> R(k, t = 0) + R(k, t = 0) (13) 

In Eq.(p"3]), R(k, t) describes the reflection coefficient corresponding to the initial condition. Since the relation between 
the reflection coefficient and the potential function is nonlinear, the perturbation of the reflection coefficient R(k, t = 0) 
cannot be associated with the spectral reflection coefficient corresponding to the initial condition fo(x). However, we 
can construct any initial condition fo(x), by imposing special conditions on the spectral reflection coefficient R(k,t). 
With this we mean, the we can compute and analyze the time-evolution of different classes of perturbations f(x,t), 
by changing the analytical structure of R(k,t). If both the unperturbed reflection coefficient and the perturbation of 
the reflection coefficient are rational functions of the wave-number, analytic expressions for f(x,t) can be obtained. 
Suppose R(k,t) is a spectral reflection coefficient that can be associated with the unperturbed solution u(x,t). In 
Appendix A, analytical expressions for the solution u(x,t) are derived if the reflection coefficient R(k,t) is a rational 
function of the wavenumber. From the unperturbed reflection coefficient R{k 1 1), we can construct a kernel K(x, x, t) 
(Appendix A): 

*<*•*•" = lH ' 14 » 

In Eq.([l4|) stands the prime for the derivative with respect to the space-coordinate x. The determinant T>{x,t) in 
Eq.(|l4|) is given by: 

V(x,t) = det {fy - (p l + Pj r 1 R J e 2 ^ x+i ^} (15) 

In Eq. (|l5|) , pi are the poles and Ri the residues of the unperturbed reflection coefficient R(k,t) Solutions of the 
KdV-equation can be derived by taking the following derivative: 

u(x, t) = -2^-K(x, x, t) (16) 
ax 

If the reflection coefficient R(k,t) is contaminated with a small perturbation R(k,t), Eq. (|l5|) contains the poles 
and residues of both the unperturbed reflection coefficient R(k, t) and the reflection coefficient corresponding to the 
perturbation R(k,t). It is shown in ref. j6| , that if the reflection coefficient R(k,t) undergoes a perturbation as given 
in Eq. (|f3|) , the determinant T)(x,t) can be expanded into the following series: 

V{x,t) -> V(x,t) + £ (x,t) (17) 

In Eq. (p~7|) , T>(x, t) is the determinant (|l5|) in absence of perturbations. The effect of the perturbation is expressed in 



the determinant £{x,t). If Eq.(17) is substituted in Eq.(|14|) we obtain the following result: 



K(x,x,t)^^4±^4 (18) 
v ' ' ; V{x,t) +S(x,t) v ; 
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Using some basic algebra, the kernel associated with the unperturbed solution u(x.t) can be separated from the 
right-hand side of Eq. : 

V(x,t) V(x,t)£'(x,t)-V'(x,t)£(x,t) 
{,,) V(x,t) + V(x,t)[V(x,t) + £(x,t)\ { ' 

The term V'(x,t)/V(x,t) in Eq.(|l|) can be identified with the time-evolution of the unperturbed problem. From 
Eq.([l^), we can identify an expression for the perturbation f(x,t): 

nd { V(x,t)£'(x,t)-V'(x,t)£(x,t) \ 
I[ ,j ~ dx\ V(x,t)[V(x,i)+£(x,t)} J [ > 

It should be realized that in the determinant £(x,t) the poles and residues of both the reflection coefficients R(k,t) 



and R{k,t) are present. It follows from Eq.(19) that the perturbation f(x,t) is large with respect to the unperturbed 
solution if the denominator T>(x,t)[V(x,t) + £(x,t)] of Eq.(pO|) is small. As illustrated in the following examples, in 
this case the perturbation /(x, t) can dominate the total solution of the KdV-equation. 

In the following examples, we consider the case in which the unperturbed solution of the KdV-equation has one 
single pole (p = i(3) and one residue (R = id). The unperturbed solution of the KdV-equation has either soliton-likc 
behavior as in Eq.(^) or singular behavior depending on the position of the pole and the residue. In the following 
experiment, we choose in case of the soliton d = — 1 and = 1. We can illustrate the effects of perturbations on 
the soliton by adding a certain number of poles and residues to the unperturbed determinant T>(x,t). In the lower 
panel of Fig. 3a, the time-evolution of a contaminated soliton u(x, i) is plotted. As a reference, in the upper panel of 
Fig. 3a the time-evolution of the unperturbed soliton is given. It can be observed from Fig. 3a that the effects of the 
contamination either spreads out, or travel with a different velocity. This implies that after a certain amount of time 
the unperturbed soliton and the effects of the perturbation separate. This is more clear in Fig. 3b. The short-dashed 
line in Fig. 3b indicates the initial condition and the solid line indicates unperturbed initial condition. The long-dashed 
line indicates the solution at t = 1. In the example of Fig. 3, the positions of the poles and residues is chosen in such a 
way that the numerical value of the denominator of Eq. ( pp| ) does not differ significantly from the numerical value of 
T>(x,t). As a result of this the perturbation remains in the same order of magnitude as the unperturbed solution. 

The situation dramatically changes if we chose d = 0.01 and 0=1 and keep the positions of all the other 9 poles 
and residues which specify the perturbation constant. By choosing d = 0.01 and = 1, the unperturbed determinant 
T>(x, t) can by equal to zero for certain values of x and t. Because of the analytical structure of the determinant (|l7|), 
the perturbation £ (x,t) has in this case also zeros for certain values of x and t. The time-evolution is in this case 
plotted in the lower panel of Fig. 4. This result has to be compared with the time-evolution of the unperturbed case 
which is given in the upper panel in Fig. 4. It follows that in this special case, the perturbed solution consists of two 
different branches. The first branch propagates with a similar propagation velocity as the unperturbed solution but it 
has large amplitude fluctuations on the characteristic of the unperturbed solution. The second part propagates with 
a different velocity than the unperturbed solution. In the following section, we derive analytical expressions for these 
different parts of the perturbation. We want to remark that a small perturbation of a singular solution at t = can 
always be constructed by choosing the positions of the poles and residues properly. It follows from the structure of 
equation (|2p| ) that these small perturbations can grow when V{x, t)[D(x, t) + £(x, t)] is small. 

From the results in this section we can conclude that the soliton exhibits stable behavior whereas the singular 
solution is unstable. Eq.(|20|) gives a general expression for the perturbation of a solution of the KdV-equation. The 
perturbation is small with respect to the unperturbed solution, if the unperturbed solution has no poles close to the 
origin in the complex plane since the dominator T>{x 1 t) + £ (a;, t) can not be zero. In contrast to this, the perturbation 
is large with respect to the unperturbed solution if the dominator T>{x, t) + £(x, t) in Eq.(po|) is nearly singular. This 
is the case if the unperturbed problem has poles close to the origin in the complex plane. This means that the singular 
solution is sensitive for noise. Moreover, for both the soliton and the singular solution the propagation velocity of 
the perturbation, which determines whether the perturbation separates from the unperturbed solution, is a crucial 
parameter for the stability. The propagation velocity and the amplitude of the perturbation depend of the position of 
the poles and residues of the perturbation. In the following section we study this in more detail. 



5 



4 A series solution for f(x,t) 



In order to analyze the behavior of the perturbation f(x,t) it is convenient to formulate this function as a series 
solution. Our starting point is the Marchenko equation without bound states in the wave-number domain as given in 
Appendix A. If the reflection coefficient undergoes a perturbation (|l3|), we find the following relation. 

1 . f°° .{R{k',t = 0) + R(k',t = 0)\ F{k , ,x,t)exp[2i(k'x + 4:{k'}H)} 



1 f°° 
F(k,x,t) = 1-1 lim / dk 



1+ / C(k,k',t)F(k',x,t)dk' + / C{k,k' ,t)F(k' ,x,t)dk' (21) 

J — OO J — OQ 

The function F(k,x,t) is related to the kernel K(x,y,t) by the following Fourier transform: 

POO 

K(x,y,t) = (27T)- 1 / dfce-^- x )(F(fc,z,i)-l) (22) 



Furthermore, the kernel C(k', k, t) in Eq. (pi]) is given by: 



C(k,k',t)= lim ^ — (23) 

v ; e^o+ 2in k' + k + ie v ' 

The kernel C(k', k,t) in Eq.@ is defined by: 

C(k,k',t)= lim ' , ; Uje , (24) 

v ' ' ' £ ^o+ 2in k' + k + ie y ' 

Eq.(|2l]) can be represented schematically using a Dyson's representation as given if Fig. 5. Using an iteration technique, 
the Dyson's series of Fig. 5, can be expanded. The result is given in Fig. 6. It is observed from Fig. 6, that the total 
expression for the function F(k, x, t) consists of three parts. The first part can be identified with all diagrams consisting 
of solid dots only. This series of diagrams represents the time-evolution of the unperturbed solution. The second series 
consists of diagrams having open dots only. This series of diagrams represents the time-evolution the perturbation in 
absence on the unperturbed solution. The remaining diagrams consist of a combination of both solid and open dots. 
This series of diagrams represents the interaction between the unperturbed solution and the perturbation. We can 
formally solve the Dyson's equation by iteration. The solution is shown by the diagrams in Fig. 6, and is given by: 



F(k,x,t) = l+ C(k',k,t)dk' + / C(k,k',t)C(k',k",t)dk'dk" 

J —oo J —oo J —OO 

poo pOO 

C(k',k,t)dk' + / / C(k,k',t)C(k' \k" ,t)dk'dk" + ■■■ 



— oo J — oo 



C(k,k',t)C(k',k",t)dk'dk" + / / C{k,k' ,t)C{k' ,k",t)dk'dk" + ■ ■ ■ (25) 



If both R{k, t = 0) and R{k, t = 0) are rational functions of the wave-number, the integrations in the series ([25]) can 
be carried out analytically by performing a contour integration in C + . This is justified by the fact that the reflection 
coefficient R(k,t = [))—► 0(l/k) if k — > oo. The poles of the denominator of Eq.(|25|) are all situated in (CP so the 



only contribution to the integrals of Eq.(25) comes from the poles of R(k, t = 0) which are situated in C . We write 
Eq.(|2|) in the following way: 

F(x, k,t)-l = F un (x, k, t) + F pe (x, k, t) + F int (x, k, t) (26) 

The function F un (x,k,t) only contains contributions of R(k,t) and can be identified with the time-evolution of the 
unperturbed solution (solid-dot diagrams). The function f(x,t) consists of contributions of both F pe (x,k,t) and 
Fint(x,k,t). The term F pe (x, k) only contains contributions of R(k,t) and can be identified with the time-evolution 
of the perturbation in the absence of u(x,t) (open-dot diagrams). The term Fi nt (x, k,t) contains contributions of 
both R(k,t) and R(k,t) and represents the interaction between the unperturbed solution u(x,t) and the perturbation 
(diagrams consisting of solid-dots and open-dots). In the following we analyze these three terms separately. 
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4.1 The time-evolution of the non-interaction terms 

Suppose the unperturbed reflection coefficient R(k,t) has N poles. It follows from Eq.(|2^) and Eg. (|26|) that the 
contribution to the function F(k,x,t) from the terms containing only the reflection coefficient i?(fc,i) (solid-dot 
diagrams) is equal to: 



N R N R R 

F un (x k t) = V n% c 2i(PiX+ipft) + n i n ] e 2i{(p f +p 3 )x+4(p?+pf)t} 



AT 

R i R 3 R l e 2i{(pi+Pj+Pi)^+4(pf+p 3 3 +pf)t} . . . (- 2 7) 



where pi and i?i are the poles and residues of the unperturbed reflection coefficient. If the Fourier transform ( J22Q is 
now performed on the unperturbed part of F m [x, k, t), we obtain the following expression for the kernel K un {x, y, t): 

N N R R 

K un (x,y,t)=iYR i e 2ip ^+v) +Sip * t + i V -i^ e ^ + ^V^e 8 ^ + ^ t + 

N 

i \^ RjRjRi e ip l (x+ v ) e 2i(p i + Vj )x e si( P ?+p 3 j +pf)t ! (28) 

ij^i(Pi+PO(Pi+!>i) 



After putting y = x and taking the derivative: 



u un (x,t) = -2-^-K un (x,x,t), (29) 
eta; 



the following expression for the unperturbed solution is obtained. 

N N 



u un (x,t) = 4^%e 2 *- 1+4 ^+4^ RiRje 2i{(Pi+Pi)*+*(Pi+Pm + 
i=l «,3=1 

4 y^ (RjRjRQiPi +Pj + pi) ^iUPi+Pi+piU+MpS+p^+pfK} ( 30 ) 

This result is already obtained in ref. Q . From this result we can conclude that a general solution of the KdV-equation 
can be expanded in an infinite series of exponential basis functions. 

In a similar manner we can derive the time-evolution of the contributions to the Dyson's series in Fig. 6 for all 
the terms that can be identified with the perturbation only (open-dot diagrams). Suppose the perturbation on the 
reflection coefficient R(k, t) has M poles, it follows using a similar argumentation as for the evaluation of the term 
u un (x, k, t) that: 

M M 

u pe (x,t) = 4^fi ! p ! e 2 * I+4 ^+4^ Rjl. e 2im+Pj)x+^+P 3 j )t} + 

i=l = l 



4 y^ (RjRjRi)(Pi +Pj +pi) e 2i{(p,+p : ,+p i )x+4(p?+^+pf)t} + . . . ( -.) | j 

ijjLl iPr+P 3 ){Pj+Pl) 

We observe from this result that both the unperturbed solution and the time-evolution of the non-interaction elements 
have a similar analytical structure. We remark that the solution u un {x,t) corresponds to a series expansion of the 
terms in Eq.([l9|) is which only the determinant T>(x, t) is present. We can conclude that both the unperturbed solution 
and the perturbation evolve as an infinite series of non-dispersive solutions in time. However, the spectral components 
of the perturbation u pe {x, t) generally travel with a different velocity than u un (x, t). This is already observed in Fig.l 
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and Fig. 2. In Fig.lc it is shown that a perturbation of a nondispersive solution of the KdV-equation introduces a 
dispersion effect, and in Fig. Id, it is shown that this results in a different propagation velocity of perturbation. If the 
unperturbed solution is a localized non-dispersive function, both the unperturbed solution and the perturbation travel 
with a different velocity depending on the position of the poles of the perturbation. The velocity of every spectral 
component of the perturbation is determined by its corresponding pole position. In the following we examine the 
behavior of the interaction terms. 

4.2 The time-evolution of the interaction terms 

The function Fi n t(x,k,t) corresponding to the interaction term contains contributions of both the unperturbed reflec- 
tion coefficient R(k, i) and the perturbation of the reflection coefficient R(k, t). In Fig. 6, it consists of all the diagrams 
having both solid and open dots. An analytic expression of all these diagrams is given by the following equation: 

/OO />00 rOO /"OO 

/ C(k,k',t)C{k',k",t)dk'dk" + / / C(k,k',t)C(k',k",t)dk'dk" 
- oo J — OO J — oo «/ — OO 

OO pOO pOQ 

/ / C(k,k',t)C(k',k",t)C(k",k"',t)dk'dk"dk"' + ■■■ (32) 

— oo J — oo J —oo 

Since, both R(k, i) and i?(fc, t) are rational functions of the wave-number, we can carry out the integrations in Eq.(|32]) 
analytically. If we proceed in a similar manner as for the non-interaction terms, we obtain the following result for the 
time-evolution of the interaction terms: 

jV M 

u int (x,t) = 8j2J2^e 2i{(pi+ ^ yx+4{p!+ ^ )t} 

i=l 3=1 

+ 12 V V RiR j Rl{ - p * tli + P l \ ^(P,+P,+P,)x+M P ?+pUpf)t} (33) 

hki (Pi+pM+pi) 



3, 

N M 



R % R 1 Rl{p i +pj 2i{(p i +p J +p i )x+4(pf+pf+pf)t} 



(P* + Pi)(Pi + Pi) 



As witnessed from the previous section, it follows from Eq.(|33[) that the interaction term is large with respect to the 
unperturbed solution, if the unperturbed problem has poles close to the origin in the complex plane. The analytical 
structure of Eq. (|30|) , Eq.(|3l|) and Eq.(j33[) enables us to formulate an expression for the amplitude behavior of the 
perturbation f{x,t). If the perturbation does not have poles close to the origin in the complex plane and if the 
perturbation is small with respect to the unperturbed problem at t = 0, it follows from Eq. (|30|) and Eq.(|3|) that the 
term u pe {x,t) remains small with respect to u un (x,t). However, as already concluded in the previous section, if the 



unperturbed solution has singularities, the structure of the interaction term ( J33[ ) introduces necessarily singularities 
in the time-evolution of the perturbed problem. The physical meaning of Ui n t(x, t) is visualized in Fig. 4. The example 
given in Fig. 4 only differs from the example given in Fig. 3, by the position of the poles and residues that characterize 
the unperturbed solution. This implies that the function u pe (x, t) in Fig. 4 does not differ from that in Fig. 3. However, 
due to the large amplitudes in Fig. 4, u pe {x,t) is small with respect to u un {x,t). In the upper panel of Fig. 4, the 
unperturbed singular solution is plotted. In the lower panel of Fig. 4, the perturbed solution is plotted. As remarked 
in Sec. 3, in this special case, the perturbation consists of two parts. One part propagates over a different characteristic 
than u un (x, t). The second part propagates over the characteristic of u un (x,t) and is responsible for large amplitude 
fluctuations on the characteristic of the unperturbed solution. It is easy to see that an interaction term as given by 
equation ( |33| ) introduces large amplitude fluctuations. If we assume that the unperturbed solution consist of one pole 
(pi = p) and one residue, than we find at the characteristic x — —4p 2 t: 

M 

Umt ( x = -4p 2 ,t) = 8^JJfl j e 85 i^- ! '^ i +Lo.t, (34) 
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From this result it follows that the interaction term fi n t(x,t) introduces fluctuations at the characteristic of the 
unperturbed solutions. This result explains the behavior of the perturbed singular solution in Fig. 4. It is observed in 
this figure that certain spectral components of the noise travel with a different velocity, and that at the characteristic 
of the unperturbed solution large amplitude fluctuations occur. This is the result of the presence of the interaction 
term Ui nt (x,t) which has the same magnitude as the unperturbed solution u un {x,t). The interaction term is large 
with respect to the unperturbed solution if <C pj for all possible Ri and Pi. If this condition is satisfied, the 
amplitude of the perturbation f(x,t) is small with respect to the amplitude of the unperturbed solution. 



5 Numerical example 

In this section the results that are obtained analytically in the previous sections are illustrated numerically in the case 
of a reflection coefficient consisting of an infinite number of poles and residues. In this section we give an numerical 
example in to illustrate the stability of the soliton. In a discrete representation the KdV-equation takes the following 
form: 

z+i _ i | A . J u n+3 ~~ 3u l n+1 + 3u l n _ 1 — u n _ a 6u l n u* n+1 — 6u* n u l n 1 

In Eq. ( [35|) , the solution of the KdV-equation at time i At and position nAx is given by u l n . In Eq. , At represents 
the time-step and Ax represents the distance between two grid-points. The discrete KdV-equation ( |35|) is solved 
numerically using the fourth-order Runge-Kutta scheme as given in ref. || . 

In the numerical example that follows, the KdV-equation is solved on a line-segment of a total length of 87r. The 
initial condition given in Fig. 7a consists of the standard soliton which is used in the previous sections (J3 = 1 and 
d = —1). The soliton is contaminated with a noise-function having a maximum amplitude of 10% of the maximum 
amplitude of the soliton. In Fig. 7b, the soliton after a propagation-time of 1.1 sec is given, it can be seen that the 
contamination has started propagating out of the soliton. This process continues, and at t = 2.2 sec, the contamination 
has virtually propagated out of the soliton. 

This experiment is a numerical confirmation of the experiments of the previous sections. It reflects the case 
that small perturbations propagate out of the unperturbed soliton so that the noise separates form the unperturbed 
solution. This is also the reason why we observe solitons in nature. In the example of Fig. 7 the spectral contents of the 
perturbation are chosen in such a way that the soliton and all the spectral components of the perturbation travel in 
such a way that the unperturbed solution and the perturbation separate during the course of time. In a real physical 
situation solitons are always contaminated with noise at time t = 0. If we observe solitons in nature, it follows from 
this paper that the spectral components of the perturbation have a small amplitude for low frequencies. As a result 
of this, the soliton and the perturbation separate during the course of time and the soliton is "born" . If the spectral 
components of the noise function are of the same magnitude than the spectral components of the soliton, the noise 
propagates with a similar velocity as the soliton, and hence we can conclude that no soliton is created. 



6 Conclusions 

In this paper the stability of the KdV-equation is investigated. In particular, attention is paid to the stability of 
localized solutions. Is Sec. 2, the effects of nonlinearity and dispersion are investigated. For non-dispersive solutions of 
the KdV-equation these effects have to be in balance. It is pointed out in Sec. 2 that in the cases of small perturbations 
the balance the dispersion and the nonlinearity is disturbed. As a result of this an additional dispersion effect is 
introduced, which generates the "force" that separates the noise from the unperturbed solution. 

In Sec. 3, inverse scattering techniques are used to formulate an analytical expression which describes the behavior 
of perturbations of solutions of the KdV-equation. From the examples of Sec. 3, it can be concluded for soliton-like 
solutions of the KdV-equation, that stability implies that noise either propagates out of the unperturbed solution, or 
that the noise spreads out so that the amplitude is reduced. Furthermore, it is observed for singular solutions that 
although the noise is propagating partly out of the unperturbed solution, large amplitude variations contaminate the 
perturbed solution. 

This behavior is examined in Sec. 4, by expanding the perturbation into a series-solution. It is concluded in 
Sec. 4 that the behavior of the perturbation f(x,t) of the KdV-equation is strongly correlated to the behavior of the 
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unperturbed solution. If the unperturbed problem has singularities the amplitude of the perturbation is in the same 
order of magnitude as the amplitude of the unperturbed problem. This result explains the large amplitude variations 
from which the perturbed singular singular problems suffers. Furthermore, in Sec. 4, a criterion for the position of the 
poles and residues of the perturbation is given to posses a stable behavior of the soliton. 

Lastly, in Sec. 5 we give a numerical example is given to illustrate the stability of the soliton. The perturbation 
used in Sec. 5 has an infinite number of poles and residues. We observe that in the numerical case the same conclusions 
can be drawn as in cases of analytical perturbations: the soliton possesses stable behavior if the noise is propagating 
out of the soliton. 
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Appendix A: The inverse problem for rational reflection coefficients 

In this appendix a brief formulation of the inverse problem for rational reflection coefficients based upon the formulation 
of Sabatier [JtJ is given. For a detailed treatment of the mathematics we refer to the book of Chadan and Sabatier ||] . 
Our starting-point is the following equation: 

1 , f°° l-T(k',t)F=(k',x,t) „, 
F±(k,x,t)-1 = — lim / ' dk' 



2ni c—o+ J_ 00 k' + k + ie 



J_ Um f°° R±(k\t)F ± (k\x,t)expl±2ik'x} ^ 
2iri eJo+ ./ ^ k! + k + ie 



In Eq.(A-l) F + (k,x) is defined for x > 0, and F^(k,x,t) for x < 0, k G (D. The function F±(k,x,t) is defined by: 

F±{k,x,t) = exp[Tikx)f±(k,x,t) (A-2) 

The Jost solutions f±{k,x,t) are those solutions of the Schrodinger equation satisfying the following boundary condi- 
tions: 



f+{k, x, t) : lim e- lkx f+(k, x,t) = 1 (A-3) 

X — >-OC 

f-(k,x,t): lim e ikx f-(k,x,t) = l (A-4) 



They satisfy the following integral equations: 

f+(k,x,t) = e lk * - £ Sm KX k - V) V(y, t)f+(k, y, t)dy (A-5) 

f-(k,x,t) = e- ik *-J X SmkiX k ~ y) V(y,t)Uk,y,t)dy (A-6) 

It is well known that the functions f±(k,x,t) and therefore also the functions F±(k,x,t) are holomorphic in C + B. 
The potential V^(x, t) has to be in the Faddeev class L\: 

(l + ]a;|)|V(a:,t)l < oo (A-7) 

The scattering coefficients R+(k, t), R-(k, t), T(k, t) are defined by the asymptotic behavior of the physical solution of 
the Schrodinger equation: 

Vi{k,x,t) ~ 4^ T{k,t)e ikx x -» +oo (A_8J 

In the case of rational reflection coefficients they take the following form JtJ : 

jmm=o> = ,. n^nn < A - io) 
r(M = 0) - nfefe^f (A - n) 

Following Sabatier Q, the degree of the polynomial -P(fc) has to be smaller than q. Further, Im /ij > except if 
/ii = 0, Im A; < 0. The transmission coefficient T(k,t) is supposed to be an irreducible fraction, and the sets M + , 
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M _ , L + , L~ contain numbers =/= 0. If the potential is real then are pure imaginary. If /ifc,Afc are not pure 

imaginary then there exists — —X* k - It can be shown that T(k) is meromorphic in C + and if there are pol es th ey 
are in Im k |§. If there are no bound states, T(k, t)F T (k, x, t) is holomorphic in (D + and the first integr al of ( A-l ) is 
zero. If T(k,t)Fzf(k,x,t) is holomorphic in C + and all the poles pi of R+(k,t) are simple, the integral ( A-l ) can be 
solved by contour integration in the upper-half plane. The result is: 



F+(k,x,t) 



1=E 



R J F + {p j ,x,t)e 2i ^ x - 4p3 J t) 
Pj + k 



(A-13) 



The time-evolution of the residues that is used in equation ( A-13 ), is given in ref.§j. Equation ( |A-13| ) can be solved by 
letting take k the values of the discrete poles p <G V . We then obtain a linear set of algebraic equations that determine 
F + (pj,x,t) for all values of j. This set o f equations can be solved making use of Cramer's rule. We obtain after 
resubstituting the result in equation ( A-13| ), using equation ( p2| ) and putting x = y: 



K + (x, x, t) 



V + {x,t) 



where: 



V+(x,t) = dei{% - ^ +p J )- 1 R 1 e 2 ^ x ~^} 
and T>' + (x,t) is the derivative of D+(x,t) with respect to x. 



(A-14) 



(A-15) 
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Captions for figures 



Figure 1: A: Unperturbed soliton (solid line) and perturbed soliton (dashed line). B: the dispersion of the unperturbed 
soliton (short-dashed), the nonlinearity of the unperturbed soliton (long-dashed line) and the sum of the dispersion 
and nonlinearity (solid line). C: Similar as for Fig. lb, but now for the perturbed soliton. D: The time-derivative of 
the perturbation. E: Finite-difference solution of the perturbed soliton at t = 2.2 sec. 

Figure 2: A: Singular solution (solid line), the dispersion (short-dashed line) and the nonlinearity (long-dashed line). 
B: The dispersion is case of a 10% amplitude error (short-dashed line). The nonlinearity in case of a 10% amplitude 
error (long-dashed line) and the sum of the nonlinearity and the dispersion (solid- line) . 

Figure 3: A: Time-evolution of a soliton in case of a perturbation consisting of nine additional poles and residues (lower 
panel) and the time-evolution of the unperturbed soliton, (upper panel). B: The initial condition (short-dashed), a the 
solution at t = 1 (long-dashed line) and the unperturbed soliton (solid- line). 

Figure 4: Time-evolution is case of a singular solution contaminated with the same perturbation as in the previous 
figure (lower panel) and the time-evolution of the unperturbed singular solution (upper panel). 

Figure 5: Dyson's representation for the integral equation ( [H|) 

Figure 6: Dyson's representation of the series expansion of Fig. 5. 

Figure 7: A: Plot of the soliton contaminated with a perturbation consisting of an infinite number of poles and 
residues. B: The initial condition, the solution at t=l.l sec and the solution at time t — 2.2 sec. 
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